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Abstract —In this paper, we propose and study the effec¬ 
tiveness of customer engagement plans that clearly specify the 
amount of intervention in customer’s load settings by the grid 
operator for peak load reduction. We suggest two different 
types of plans, including Constant Deviation Plans (CDPs) and 
Proportional Deviation Plans (PDPs). We define an adjustable 
reference temperature for both CDPs and PDPs to limit the 
output temperature of each thermostat load and to control the 
number of devices eligible to participate in Demand Response 
Program (DRP). We model thermostat loads as power throttling 
devices and design algorithms to evaluate the impact of power 
throttling states and plan parameters on peak load reduction. 
Based on the simulation results, we recommend PDPs to the 
customers of a residential community with variable thermostat 
set point preferences, while CDPs are suitable for customers 
with similar thermostat set point preferences. If thermostat loads 
have multiple power throttling states, customer engagement plans 
with less temperature deviations from thermostat set points are 
recommended. Contrary to classical ON/OFF control, higher 
temperature deviations are required to achieve similar amount 
of peak load reduction. Several other interesting tradeoffs and 
useful guidelines for designing mutually beneficial incentives for 
both the grid operator and customers can also be identified. 

Index Terms —Smart grid, user inconvenience, peak load, 
customer engagement plan, demand response. 


Nomenclature 

K Number of operable states of thermostat loads. 

K 3 Number of operable states of thermostat load i of 

customer j. 

k Index of each operable state of thermostat load. 

J Number of customers in the residential community. 

I Set of flexible loads for which customer engage¬ 

ment plans are defined. 

It Set of thermostat loads for which customer engage¬ 
ment plans are defined. 

Is Set of shiftable loads for which customer engage¬ 
ment plans are defined. 

Itc Set of thermostat loads used for cooling for which 
customer engagement plans are defined. 

Ith Set of thermostat loads used for heating for which 
customer engagement plans are defined. 
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P Set of flexible loads owned by customer j. 

I 3 t Set of thermostat loads owned by customer j. 

I J S Set of shiftable loads owned by customer j. 

I J Tc Set of thermostat loads used for cooling owned by 

customer j. 

P Th Set of thermostat loads used for heating owned by 
customer j. 

i Each flexible load index. 

j Each customer index. 

n Index of each local peak in Algorithm 2. 

Oj^ Preference of customer j for the set point of 

thermostat load i. 

Of f Reference temperature for thermostat load i. 

A6™ x Constant value representing the maximum temper¬ 

ature deviation for thermostat load i. 

AOj^i Inconvenience severity experienced by customer j 
for thermostat load i. 

Oj,i(t) Output temperature of thermostat load i of cus¬ 
tomer j. 

Oj AC (t) Output temperature of AC of customer j , which is 
operated in state k , at time t. 

6j wh(^) Output temperature of WH of customer j , which 
is operated in state k , at time t. 

Oj : i(t) Output temperature of thermostat load i of cus¬ 
tomer j at time t during Algorithm 2 computations. 

^ac Average temperature deviation of AC loads over 
the demanded intervals in the simulations. 
Average temperature deviation of WH loads over 
the demanded intervals in the simulations. 

#iniet Temperature of the inlet water of WH. 

O a (t) Room temperature at t. 

fr(t) Rate of water flow in WH. 

Vtank Volume of the water tank of WH. 

A tank Tank surface area of WH. 

i?tank Heat resistance of the water tank of WH. 

Gj ( t ) Heat gain rate of the house of customer j at time 

t. 

T Number of intervals in the considered time dura¬ 

tion. 

t Index of each time interval. 

At Time duration of each interval. 

As™ ax Maximum scheduling delay of each shiftable load 
i. 


Ad™ ax Maximum inconvenience duration of each thermo¬ 
stat load i. 

t^tart Actual start time of shiftable load i of customer j. 
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£^ f Preferred start time of shiftable load i of customer 
3- 

tf ir Maximum number of time slots a thermostat load 

i can deny operation at its full rated power. 
t“f' n Time index corresponding to n-th largest local 
peak in the demanded operation interval of ther¬ 
mostat load i of customer j in Algorithm 2. 
fpea k Time index vector containing local peaks in Algo¬ 
rithm 2. 

hj' /L Set of time indexes during which customer j 
demands thermostat load i. 

Pi Scaling factor for thermostat load i. 
wDemand status vector of flexible device i of cus¬ 
tomer j. 

Entry, which is a binary variable, of the demand 
status vector at t. 

Pjj Rated power of flexible load i of customer j. 

I Index to denote scheduling delay in terms of 

number of time slots in Algorithm [TJ 
p j : i Vector of demanded power requirement by flexible 

load i of customer j. 

Pj,i(t) Demanded power consumption requirement of 
load i of customer j at time t. 
p l J Z Demanded power requirement vector of shiftable 
load i of customer j when delayed by l time slots. 
C jj A matrix of binary variables of order T x Kl 

Cjj(t) Binary variable indicating the operational state of 
thermostat load i of customer j in state k at time 
t. 

cj,z(£) £-th row vector of matrix C^. 

Cjj k- th column vector of matrix C^. 
a Energy required for a unit degree rise in room 

temperature. 

1(.) (•) x 1 vector of all l’s. 

e r K j Vector of all l’s except a 0 at K \-th position, 
x A given aggregated load profile vector of residen¬ 

tial community. 

x(t) Aggregated load demand at time t. 
x Aggregated load profile obtained by solving sub¬ 

problem 1. 

Xj Aggregated load profile obtained from the cus¬ 
tomer j in Algorithm Q] 
x Final aggregated output load profile. 

2F{’) A circular shift operator. 

m Particular order of shiftable devices, 

m* Optimal order of shiftable devices. 
y l Aggregated load profile in Algorithm Q] when 
shiftable load i of customer j is delayed by l time 
slots. 

pi Peak value of y l in Algorithm Q] 
l* Optimal scheduling delay for shiftable load i de¬ 

termined by Algorithm Q] 

y Output of Algorithm Q] (aggregated load profile). 

y m Output of Algorithm \T\ for ra-th particular order 
(aggregated load profile). 

pm Peak value of y m (for shiftable device order m). 

Power consumption of thermostat load i of cus¬ 


tomer j while operating in state k at time t. 

Q A matrix of order T x Kj, which has the entry of 

Vj,k,i e It. 

Zj AC Cooling capacity of customer j’s AC operating in 
state k at time t. 

EER Energy efficiency ratio of the AC. 

h Output of Algorithm 2 (aggregated load profile). 

Nac Number of eligible AC loads in the simulations. 

Nac Number of eligible WH loads in the simulations. 

I. Introduction 

Demand Response Program (DRP) can be used to reduce 
cost and improve efficiency of power grids by engaging 
customers and modifying their power consumption pattern m- 
Q. Current smart metering and bi-directional communication 
technologies, allow the inclusion of industrial as well as resi¬ 
dential consumers in DRPs lf8l- fl3l . Customer engagement for 
peak load reduction can be achieved through motivating cus¬ 
tomers by either time-based DRP or incentive-based DRP [ 14], 
aa. In time-based DRP, electricity prices are dynamically 
varied at different times and the customers can modify their 
electricity consumption in response to the change of price. 
In incentive-based DRP, fixed or time-varying payments are 
offered to the customers under specific constraints. In general, 
residential customers are not very responsive to DRPs due 
to the uncertainty in their electricity bills and the lack of 
clarity about the resulting inconvenience that they might expe¬ 
rience 03- The participation of residential customers in DRP, 
however, is extremely important to address the predictable 
and non-predictable supply and demand variations in order 
to reduce the electricity generation cost 03 , 03 . 

The increased penetration and usage of Air Conditioner 
(AC) and Water Heater (WH) loads in residential sector are the 
main reasons for increasing the peak load demand on the grid 
ED-ED. in d, the authors propose a direct load control 
algorithm to investigate the potential of WH loads in providing 
load-balancing service, while in 1231 . the authors evaluate 
the potential of Heating, Ventilation and Air-Conditioning 
(HVAC) loads for providing regulation services. These two 
papers suggest that the aggregated regulation service provided 
by the WH and HVAC loads can become one of the major 
sources of revenue for the grid. The authors in ll24l propose 
a controller for thermostatically controlled loads, which can 
be used for peak shaving and load shifting by managing the 
aggregated HVAC load shapes. The authors in lf25ll propose 
HVAC control algorithm to regulate the indoor air temperature 
inside a defined dead band by using an ON/OFF power control. 
However, the algorithm does not consider the retail price of 
electricity in the HVAC control, and therefore may fail to 
curtail load during peak price period. 

In [26], a price responsive control strategy is proposed 
for HVACs to reduce peak load. The controller changes the 
thermostat set point of HVAC loads depending on the elec¬ 
tricity retail price published in every 15 minutes. The authors 
demonstrate a significant peak load reduction with a modest 
variation in thermal comfort. In ['27’], the authors consider the 
problem of managing large number thermostat appliances with 
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ON/OFF power control for peak load reduction. They discuss 
a centralized, model predictive approach and a distributed 
structure with a randomized dispatch strategy. 

Several authors have also designed algorithms that may 
help the grid operator to reduce peak demand by controlling 
and scheduling common household appliances. For instance, 
[28] proposes an optimal and automatic residential load com¬ 
mitment framework, in which the households achieve the 
minimum payment by responding to the time-varying prices 
offered by the grid. Similarly, a Stackelberg game is studied 
in |[29l to maximize both the revenue of the grid operator and 
the payoff to each user by determining the optimal electricity 
price and the optimal price consumption. In (30|, the authors 
use a water-filling based scheduling algorithm, which does 
not require any communication between scheduling loads, to 
obtain a flat demand curve. The authors also explore the 
possible errors in demand forecast and potential incentives 
for customer participation. A genetic algorithm based energy 
resource scheduling technique, which includes day-ahead, 
hour-ahead, and five minutes ahead scheduling is proposed for 
smart grid in ED. Finally, a heuristic method called Signaled 
Particle Swarm Optimization is proposed in [32) for distributed 
energy resource scheduling in smart grids. 

These DRP algorithms certainly present feasible solutions 
to address the problem of increased peak load on the grid. 
However, designing schemes, which can attract the interest 
of significantly large number of customers to participate in 
DRP remains a major challenge ltl6l - lH8l . The reluctance 
to participate by the residential customers is due to the lack 
of clarity in specifying key information such as the number 
of times customers would be called upon to participate, the 
range of temperature variation of AC and WH loads, and 
the financial benefits associated with such participation. In 
recent years, companies like Idaho Power have introduced 
direct load control methods, such as AC Cool Credit program, 
to switch OFF the AC loads of their residential customers 
[33], where the duration of OFF time is set according to a 
pre-determined agreement. Such agreement has been reported 
to be very effective in considerably reducing the peak load. 
However, the resulting room temperature deviations during the 
OFF intervals are not specified in the agreement, which can 
result in excessively high inconvenience for some customers. 

In this paper, we propose customer engagement plans that 
clearly specifies all the key inconvenience parameters for 
shiftable and thermostat loads, and thus may encourage the 
customers to pick plans according to their behavioral and 
financial requirement^. For each shiftable load, maximum 
scheduling delay is specified in the plan. For each thermostat 
load, maximum temperature deviation from the thermostat set 
point and the maximum time duration during which the actual 
temperature deviates from the thermostat set point are clearly 
definec@. It is important to note that this work is an extension 

designing actual financial rewards or incentives for customers is out of 
scope of this paper. 

2 Customers generally cannot understand the amount of intervention by 
specifying the energy consumption limits on their loads or by announcing 
the dynamic price on energy consumption. It is much easier for them to 
comprehend the inconvenience if specified in terms of scheduling delays and 
temperature deviations of loads. 


of our previous work in j34l . In f34l . we proposed demand 
response management plans for AC loads. We modeled the AC 
load as a power throttling thermostat device and specified the 
thermostat set point, temperature deviation and inconvenience 
duration for the customers. We showed that a power throttling 
thermostat device can operate in K > 2 states, where K is 
the number of possible power states. Note that in a 2-state 
model, thermostat loads can only be turned ON/OFF for DRP, 
whereby they can operate in K different states for K > 2. For 
example, in a 3-state model, the thermostat load can be turned 
OFF, operated at 50% of the rated power or operated at full 
rated power. In |[34l . we determined the effectiveness of such 
plans and studied the impact of temperature deviation, time 
duration of inconvenience, and the impact of increasing the 
power throttling states on peak load reduction. However, the 
inclusion of thermostat set point as a plan parameter in [34], 
requires all the customers to adjust their thermostat set point 
to the same specified value and thus experience additional 
temperature deviation for the time duration as laid out in the 
plan. Such customer engagement plans, therefore, may have 
practical limitations and fairness concerns. 

To this end, we extend the work in El in this paper 
to include multiple shiftable and thermostat loads. In this 
context, we propose two new types of customer engagement 
plans: the Constant Deviation Plan (CDP) and Proportional 
Deviation Plan (PDP), which allow customers to have different 
thermostat set points. We model thermostat loads as power 
throttling devices such as in El- However, in contrast to 
El, the model adopted in this paper is quite generic and can 
allow different customers to have different number of loads 
including thermostat loads with different number of power 
throttling states. We propose algorithms that enable us to iden¬ 
tify mutually beneficial plans for both the grid operator and 
the customers with power throttling thermostats and shiftable 
loads through the introduced customer engagement plans, i.e., 
CDPs and PDPs. Through MATLAB simulations, we study the 
impact of the number of power throttling states, magnitude of 
temperature deviations, and scheduling delays on the peak load 
reduction. The impact of scheduling delays in the presence 
of thermostat loads with multiple power throttling states is 
investigated, and certain observations and comparisons of CDP 
and PDP type plans are made. The impact of CDP and PDP 
parameters in controlling and determining the number of eli¬ 
gible thermostat loads and average temperature deviations are 
also investigated. Some useful design guidelines for CDP and 
PDP plans are provided, which could be helpful in designing 
the financial incentives for participating customers. Further, 
several new and interesting research directions are identified. 

The rest of the paper is organized as follows; the customer 
load model, residential community model and customer en¬ 
gagement plans are explained in Section II, the optimization 
problem for peak load reduction and the algorithmic solutions 
are presented in Section III, simulation results are given in 
Section IV, while the paper is concluded and future research 
directions are identified in Section V. 
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II. Customer Load, Residential Community 
Models and Customer Engagement Plans 

A. Customer Load Model 

In this paper, we group the appliances, also called loads or 
devices, in a home into two categories: essential appliances, 
and flexible appliances. The power required for all the entire 
essential appliances serves as a base load, which always has to 
be provided by the grid operator. We further classify flexible 
loads into shiftable and thermostat load^]. We assume that 
the shiftable loads, such as Dishwashers (DWs) or Clothes 
Dryer (CD), require a constant power draw for a specified 
duration. We further assume that neither the power draw nor 
the operational duration of shiftable loads can be reduced. 
The customer has a preferred operation interval for each 
shiftable load that best suits his/her behavioral requirements. 
Scheduling delays in the operation of shiftable loads, can cause 
inconvenience to customers. 

Thermostat loads such as AC and WH, on the other hand, 
have to maintain an output temperature that matches the 
customer-defined thermostat set point, and any deviation from 
these set points cause inconvenience. Temperature deviation 
can be controlled in two different ways: 1) by re-adjusting 
the thermostat set point; or 2) by re-adjusting the power 
consumption of the device. During the re-adjustment of set 
point, a smart thermostat automatically adjusts the thermostat 
settings to control temperature deviations. One such example 
is Google Nest, which can be used to re-adjust the thermostat 
set point of AC loads f35l . Conversely, power consumption 
can be re-adjusted by providing a smart electrical interface 
with each device which enables it to respond to the control 
signals received from a controller. For example, Idaho Power 
controls the temperature deviations of AC loads by adjusting 
the number of ON/OFF switching frequency. Thermostat re¬ 
adjustment method requires mathematical models (for each 
thermostat load) to compute the resulting power consumption 
for some specified temperature deviation, while power re¬ 
adjustment method requires models to compute the resulting 
temperature deviations for specified power consumption. 

In this paper, we consider thermostat loads that are capable 
of operating in K > 2 power consumption states. We have 
also adopted the power re-adjustment method to control the 
temperature deviations of thermostat appliances m-m, 
l33l . |36|. This method allows us to determine the exact 
amount of power consumption by the appliances and using the 
thermal models readily available in the literature lf36l , we can 
determine the output temperature (e.g., room temperature or 
hot water temperature). Moreover, the thermal models can be 
easily extended to include multiple power throttling states. The 
AS/NZS 4755 standard (jointly adopted by the Australian and 
New Zealand Governments) ED for thermostat loads being 
manufactured and sold in Australia and New Zealand also 
mandates physical/electrical interface as well as mandatory 

3 Other categories of flexible loads e.g., electrical vehicles, which require an 
unfixed amount of power and operation duration depending on the customer 
usage are also possible. Such load categories are not considered in this paper. 
However, our framework can be easily extended in future to accommodate 
electrical vehicles or any other new category of flexible loads. 


and optional modes, which permit these loads to operate in 
different power consumption states. Variable frequency drives 
and variable speed drives are generally used to throttle power 
of thermostat loads between different states. 

B. Residential Community Model 

We model a residential community comprising of J cus¬ 
tomers (also called homes, users or consumers). We consider 
the aggregated power consumption profile of the community 
for a duration of 24 hours, and divide the total duration into 
T equal intervals. Each interval At comprises of At = 24 * 6Q 
minutes. Each customer has a base load, a set of thermostat 
loads and a set of shiftable loads. We further differentiate 
thermostat loads into two types depending on whether the 
load is used for cooling application or heating application. 
For example, in summer AC is used to cool the room, while 
WH is used to heat water. In our paper, we assume that 
the grid operator identifies a set of flexible devices denoted 
by /, further categorized into Is shiftable and It thermostat 
loads, which it wants to control. Let It c and Ith respectively 
denote the set of thermostat loads used for cooling and heating 
applications (It = It c U Ith )• The grid operator will then 
propose customer engagement plans only for these flexible 
devices. Let P C /, I° s C I s and I J T C I T respectively 
denote the set of flexible, shiftable and thermostat loads of 
customer j. Similarly, let I J T = I J Tc U I° Th , where, I J Tc C I Tc 
and P Th C I T h denotes the set of thermostat loads of customer 
j used respectively for cooling and heating applications. All 
the remaining devices of customer j are treated as essential 
loads. For example, if the grid operator is only interested in 
controlling the AC load in a residential community and defines 
a customer engagement plan only for this load, then all the 
remaining loads such as CD or WH of all the customers will 
be treated as essential loads. This model, however, is flexible 
enough to allow different customers to have different number 
of flexible devices. This model also allows P to be a strict 
subset of /, i.e., it is possible for some customers to own less 
flexible loads than the number of flexible devices / for which 
the grid operator has defined customer engagement plans. In 
such a case, it is up to the grid operator to either allow such 
customers to participate in DRP with their available devices 
and offer them reduced financial incentives or does not allow 
them at all. In the rest of the paper, in order to have a generic 
model, we assume that a customer with P c I is allowed by 
the grid operator to participate in the DRP (the problem on 
how the financial incentives are computed is out of scope of 
this paper). 

Each home is assumed to be equipped with a home con¬ 
troller as part of its advance metering infrastructure, which acts 
as an interface between the customer and the grid operator. 
There is a two-way communication link between the grid 
operator and each home controller. In our proposed method, 
the grid operator and the home controller only exchange the 
aggregated load profiles in order to minimize the commu¬ 
nication bandwidth required and also to protect the privacy 
of individual customers. The home controller executes the 
algorithms (proposed in Sections III-A and III-B) in order 
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to determine the scheduling time and power throttling states 
of its appliances according to the customer engagement plan. 
These values are stored at the home controller and the new 
aggregated load profile is communicated back to the grid. The 
home controller then schedules and operates its appliances 
accordingly. There is a two-way communication link between 
the home controller and every flexible load (shiftable or 
thermostat). The power flow is unidirectional i.e., from the 
grid operator to the home controller and appliances. The 
home controller generates appropriate control signals for each 
of its flexible device. These devices report back their status 
(ON/OFF status, power throttling state, temperature etc.) to 
the home controller. The implementation and communication 
requirements of our proposed methodology are further elabo¬ 
rated in Fig. Q] 



* 


Fig. 1. System Model elaborating the implementation and communication 
requirements of our proposed methodology 


C. Customer Engagement Plans 

Customer engagement plans can be designed in one of the 
following ways: 

• by defining the power consumption pattern for customers 
e.g., by announcing power consumption limits in peak 
load hours; 

• by defining the inconvenience parameters for flexible 
loads of the customers. 

Generally, it is difficult for customers to translate power 
consumption limits into scheduling delays and temperature 
deviations experienced by the loads, and easier for the grid 
operator to understand the amount of peak reduction and 
demand shaping. On the other hand, customer engagement 
plans that define the inconvenience parameters for the flexible 
loads can be easily understood by the customers, while their 
effectiveness for the grid operator is unclear. The second 
approach, however, is more promising to convince and en¬ 
courage customer participation (because it provides clarity to 
the customers), which is therefore adopted in this paper. 

In our proposed customer engagement plans, the inconve¬ 
nience related to any shiftable load i G Is is defined in terms 
of scheduling delay from the preferred start time of the load. 


The plan defines As™ ax as the maximum scheduling delay 
(in minutes) for each shiftable device i. The inconvenience 
related to any thermostat load i G It have two dimensions: 
1) inconvenience duration: which is defined as the total time 
duration during the demanded interval the thermostat load is 
denied operation at its full rated power; and 2) inconvenience 
severity: which is defined as the temperature deviation from 
the desired thermostat set point. Maximum inconvenience 
duration Ad™ ax (in minutes) for each thermostat load i is 
specified in the plan. However, defining inconvenience severity 
is not straightforward since each customer has a different pref¬ 
erence for thermostat set point (denoted by Oj ^), which also 
impacts the experienced inconvenience and thermal comfort 
ED, l39l . For example, if the outdoor temperature on a hot 
summer day is 82°F, then a temperature deviation of 4°F is 
more severe for a customer with AC thermostat set point at 
76°F compared to another customer with AC thermostat set 
point at 70°F. We therefore propose two different types of 
plans: CDPs and PDPs. For each thermostat load z, both CDP 
and PDP define a reference temperature denoted by 0 f f . 

In CDPs, a constant value representing the maximum tem¬ 
perature deviation for each thermostat load i, denoted by 
A 0™ x , is announced. Accordingly, the inconvenience severity 
experienced by customer j for thermostat load i is computed 
as, 

[ (mm( 6 f -0 jti ,A6^)) + V* e it. 

= < ) ; <+ (i) 

{ (mm(0 jti -0F,A0?**)) Vi e If 

where, x + is defined by, x + = max{0, x}. In PDPs, a scaling 
factor 0 < ^ < 1 is defined in the plan for each thermostat 
load i . The inconvenience severity experienced by customer j 
for thermostat load i is then computed as, 


A 9 




A x (of - + 

off 


Vi e If 

^ ^ 1Th 


( 2 ) 


Pi x (Oj, 

To summarize; for each thermostat load i following parameters 
are defined by the grid operator: 

• CDP: (Ad™ x ,A6™ X ,6f f ) i.e., (Maximum inconvenience 
duration, Maximum temperature deviation, Reference 
temperature) and CD; 

• PDP: (A d™ x ,/3i,6f f ) i.e., (Maximum inconvenience du¬ 
ration, Scaling factor, Reference temperature) and ©. 

As mentioned before, the reference temperatures in CDPs and 
PDPs bound the output temperature of thermostat loads, and 
can also be used to control the total number of loads experienc¬ 
ing temperature deviations. For example, if the thermostat set 
point of AC is below the reference temperature, then the output 
temperature for this customer cannot exceed the reference 
temperature. On the other hand, if the thermostat set point of 
AC is above the reference temperature, then such a customer 
will be allowed to operate at the preferred set point without any 
intervention. Generally, a high value of reference temperature 
for AC and a low value of reference temperature for WH will 
make a high number of AC and WH loads eligible for DRP. 
The values of reference temperatures can be decided by the 
grid operator depending on the local weather, geographical 
location, user preferences and the desired peak load reduction. 
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A sample CDP and PDP for a residential community 
comprising of AC, WH, CD and DW as flexible loads is given 
in Table U In this table, we also give an example of resulting 
temperature deviations of AC and WH loads of two customers 
with different thermostat set points to explain the difference 
between CDP and PDP and the role of reference temperatures 
(temperature deviations are computed using © for CDP and 
© for PDP). Note that the thermostat set points of AC and 
WH loads of the customers in this example, are assumed to 
lie in a typical range associated with the thermal comfort of 
humans (38), |[39l . 

The customers may subscribe to a plan on a daily, weekly or 
monthly basis. Similarly the grid operator can also change the 
set of its customer engagement plans on monthly, bi-monthly 
or seasonal basis (depending on the geography, weather pat¬ 
terns etc.). Such issues, however, are not discussed in this 
paper and are left as an interesting future work. 

III. Optimization Problem Formulation and 
Algorithm Development 

The objective of this paper is to determine the effectiveness 
of a given customer engagement plan (CDP or PDP) for peak 
load reduction. In this section, we mathematically formulate 
the optimization problem. Let w^ = • • • 5 w j,i(T)] r 

denote the demand status vector of flexible device i of user 
j ([-] T denotes the transpose operation). Each entry 
of this vector is a binary variable: a ‘V indicates that the 
device operation is demanded, while a ‘O’ indicates that it 
is not demanded. Let pj • denote the rated power of flexible 
load i of customer j. Then p^ = [p^(l),... ,p^(T)] r = 
Pj denotes the demanded power requirement vector of 
flexible load i of customer j , while Pj,i(t ) is used to denote 
the demanded power consumption requirement of load i of 
customer j at time t. The operation of shiftable loads can 
be delayed from the preferred time interval. Let denote 
the actual start time of shiftable load i of customer j, then 
according to the plan, we have the following constraint, 

tf« - /^ r < A.s"' ax . Vj,iel 3 s (3) 

where, denotes the preferred start time of shiftable load 
i of user j. This constraint is sufficient to describe the terms 
of engagement as laid out in the plan for the shiftable loads. 
Also note that this constraint only allows shiftable loads to 
be delayed from their preferred start time (our framework 
can easily allow advance scheduling of shiftable loads, which 
however, is not considered in this paper). 

We assume that thermostat load i of customer j has the 
capability to operate in K\ possible states. This model allows 
different customers to own thermostat loads with different 
number of operable states. For example, the AC of customer 
1 can operate in only 2 states, while that of customer 2 can 
have 5 states. Let C^ denote a T x K? matrix of binary 
variables. The entries of this matrix are given by variables 
Cj-(f), which represent the operational state of thermostat load 
i G I° T of customer j in state k at time t. If thermostat load of 
customer j is operated in state k in time interval £, then c^(f) 
is set to ‘1’; otherwise to ‘O’. In the subsequent discussion, 


4 = [^(1),... ,4(T)] r will be used to denote the k- th 

.fC 

column vector, while c= [cC(f) ? ..., c- \ (t)] T will be 
used to denote the t-th row vector of the matrix C^. Since a 
thermostat load can operate in only one state in any given time 
interval t , we have the following constraint on its operation, 

C jti l K s=l t, Vj,ie4 (4) 

where, 1 K j is a K\ x 1 vector of all l’s, while 1 t denotes a 
T x 1 vector of all l’s. Similarly, we can represent the con¬ 
straint on the maximum inconvenience duration, i.e., Ad™ x 
of thermostat load i mathematically as, 

e L < A df ax , Vj, i e I° T (5) 

Y i 

where, = [1 1...1 0] i.e., it is a vector of all l’s 

except a 0 in the K\ -th position. This constraint bounds the 
number of time slots in the demanded interval during which 
the operation of thermostat load is denied at its full rated 
power. The output temperature of a thermostat load at any time 
t depends on its operational state. Let denote the 

output temperature of thermostat load i of customer j at time 
t (appropriate mathematical models are required to determine 
the output temperature as a function of power state of the load). 
During the time slots when the thermostat load is demanded 
by the customer, the difference in output temperature between 
any two consecutive intervals should be less than or equal to 
A Qjj (computed using dTJ for CDP and © for PDP), which 
can be expressed in terms of the following constraint, 


) 0j,i{c j: i(t)) - 6j,i < A 6j t i, Vj, £, i e I J T (6) 


As stated before, our objective is to determine the ef¬ 
fectiveness of customer engagement plan in terms of peak 
load reduction. Let x = [x(l)...x(T)] r denote the given 
aggregated load profile vector of the residential community, 
where, x(t ) = Then we have the following 

optimization problem, 


min max x (7) 

Co i,t start t 

0, % i 

subject to constraints: ©, ©, © and © 

The objective function © seeks to minimize the maximum 
value of x. The optimal solution of this problem depends on 
the order in which customers are considered for optimization. 
Due to the combinatorial nature of the problem, it is NP-hard 
and therefore finding the optimal solution has exponential time 
complexity m-m. We develop a sub-optimal heuristic by 
decomposing the problem into two sub-problems based on the 
load categories, i.e., shiftable and thermostat loads as shown 
in Fig. [21 We first minimize the peak load by re-scheduling the 
shiftable loads subject to constraint ©. The new aggregated 
load profile is denoted by x. The peak of this load profile is 
further reduced by controlling the thermostat loads subject to 
constraints ©, © and © to obtain the final aggregated output 
load profile x. It should be noted that the order in which these 
steps are carried out can be interchanged but could result in a 
different amount of peak load reduction. 
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TABLE I 

Sample Customer Engagement Plans and an example of Resulting Temperature Deviations in CDP and PDP for two customers 



AC 

WH 

CD 

DW 

Customer-1: AC 

Customer-1: WH 

Customer-2: AC 

Customer-2: WH 

CDP 

Ad™ 

A6 ac 




aref 

A S ™ 


01, AC 

A0i ; ac 

02, AC 

A02,AC 

01, WH 

A#i,wh 

02, WH 

A02,WH 

30 min 

4°F 

80°F 

30 min 

8°F 

108°F 

60 min 

60 min 

70°F 

4°F 

76°F 

4°F 

114°F 

6°F 

104°F 

0°F 

PDP 


5 AC 

/Q ref 

U AC 


Avh 

0 ref 

As™ 

a „max 
aS DW 

#1, AC 

A0i,ac 

02, AC 

A02,ac 

01, WH 

A#i,wh 

02, WH 

A02,WH 

30 min 

0.6 

80°F 

30 min 

0.8 

108°F 

60 min 

60 min 

70°F 

6°F 

76°F 

2.4°F 

114°F 

4.8°F 

104°F 

0°F 


Our sub-division of the problem and the subsequent algo¬ 
rithms based on the load categories provides flexibility and 
easily allows us to include or exclude some load category from 
optimization. For example, the grid operator can implement 
the framework without the shiftable load category or without 
the thermostat load category if required controls (interfaces) 
are not available. Moreover, if a third category of flexible 
loads (e.g., electrical vehicles) is available, then the framework 
can be easily extended by adding a third algorithm for the 
new category of flexible loads. Such flexibility can also be 
helpful in the actual implementation, where the proposed 
framework can be implemented in stages, as the interfaces 
and controls can be different for different load categories. 
Other alternatives, in which all the appliances of a customer 
are considered before moving on to the second customer are 
less flexible compared to our approach. 


Sub-Problem 1 Sub-Problem 2 



min max x 

v 

min 

max X 

X _> 

t sta. r t t 

.A. 

Cj i 

t 


subject to: (3) 


subject to: 

(4), (5) & (6) 


Fig. 2. Decomposition of the Optimization Problem according to Load 
Categories 

In Sections III-A and III-B, we develop distributed of¬ 
fline algorithms for the two sub-problems. We assume prior 
knowledge on the aggregated load profile of the residential 
community. This information can be obtained either from 
past power consumption patterns of the community during 
the same time period or it can be deduced using some load 
prediction/forecast models sa, in. 

Remark: It should be noted that the proposed algorithms 
are heuristics since there is no optimization on the customer 
order. One way to determine the optimal customer order is 
by exhausting all the possible options, which unfortunately is 
not possible in polynomial time for a residential community 
comprising of several hundred customers (due to NP-hardness 
of the problem). Thus, for any given customer order (among 
the huge number of possibilities), our algorithms and frame¬ 
work could be used to study the effectiveness of customer 
engagement plans and to obtain useful guidelines for designing 
these plans. 

A. Sub-Problem 1: For Shiftable Loads 

We develop a distributed algorithm in which the grid opera¬ 
tor provides coordination among the customers. To be mindful 
of the security and privacy concerns of the customers, only 
the aggregated load profiles are exchanged between the grid 


operator and customers. The given aggregated load profile of 
the community, i.e., x, is communicated to the home controller 
of customer 1, which runs an algorithm developed for shiftable 
loads (which will be explained below). This algorithm de¬ 
termines the scheduling time slots for the shiftable loads of 
customer 1 with the objective of peak load minimization and 
constraint The new aggregated load profile denoted by xi 
is sent back to the grid. The grid operator communicates xi 
to the home controller of customer 2. This process is repeated 
and the aggregated profile obtained from the last customer 
J in the sequence is denoted by xj, which is also the final 
aggregated output profile x for sub-problem 2. 

1) Algorithm for Shiftable Loads: We now explain the 
algorithm for shiftable loads that is carried out at the home 
controller. The starting time of the shiftable load i of customer 
j can lie anywhere in the interval G tjf + 

1,..., + As” 1 ”*}. Thus, when a home controller j receives 

an aggregated load profile Xy_i from the grid, the algorithm 
determines the best starting time for each of its shiftable load 
in order to minimize the peak load. Let us define a circular 
shift operator denoted by which re-arranges the entries of 
a vector by shifting them one unit to the right and moving the 
last entry to the first positiorQ. The scheduling algorithm for a 
specified sequence of shiftable devices at the home controller 
of customer j is given as Algorithm |T] 


Algorithm 1 Algorithm For Shiftable Loads 


Initialize: pL = p,Vz e I J S , y = xj -1 - £, 


for i G I° s do 

for 1 = 1: As” 

y l = y + pL 


ei 3 s P 


do 


Find the peak: fi = max y l 

p^ 1 = Ap j,i) 

end for 

Find the index: l* = min^ : 


t 


start _ j-Pref 


h: • • • , s™*} 


= + 7* — 1 

j,i fj,i ' L i 1 

y = y 1 * 1 


end for 


The optimal device order, which shiftable loads should be 
considered for scheduling, can be determined by exhausting all 
the possibilities, and this requires Algorithm []] to be executed 
(I 3 S )\ times. Let the index m denote a particular order of 
shiftable devices. For this order we denote the final output of 

4 We assume that the demanded scheduling parameters of all the shiftable 
loads are defined in such a way that despite being delayed they can always 
be completed until midnight i.e., we do not allow tasks to spill over to next 
day. 
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Algorithm \T\ by y m . Let </> m = max y m denote the peak load. 
Then the optimal device order denoted by m* is the one that 
corresponds to the minimum value in the set {0i,..., 
and the aggregated output profile Xj = y m * is communicated 
back to the grid. Typically the number of shiftable devices in 
each home is small (1 to 5 devices), which makes it possible 
to determine the optimal device order and hence the optimal 
output profile. 


B. Sub-Problem 2: For Thermostat Loads 

For thermostat loads, we again develop a distributed se¬ 
quential algorithm coordinated by the grid. The grid operator 
communicates the aggregated load profile x that is obtained 
from the solution of sub-problem 1 to the home controller of 
customer 1. At the home controller we have an algorithm for 
thermostat loads (which will be explained below) and the load 
profile obtained from this algorithm denoted by xi is sent back 
to the grid operator for onward transmission to customer 2 and 
the process is repeated until we obtain the final aggregated 
output profile x. 

The power consumption of thermostat load i of customer 
j at time t when it is operating in state k can be modeled 
in terms of the rated power of the appliance by the following 
equation, 

<&(*) = fff x Pli, Vi, Me 4 ( 8 ) 

Kj -1 

For example, an AC that can operate in 5-states has the 
capability to throttle power at 0%, 25%, 50%, 75% and 100% 
of the rated power. Let us define a T x K° { matrix Q^, where 
the entries gj^(t) of this matrix are given according to ®. 

We also require thermal load models that can relate the 
output temperature obtained by an appliance when it is op¬ 
erated in some power throttling state. Each thermostat load 
has its own thermal characteristics and hence the modeling 
varies. Since AC and WH are the most significant thermostat 
loads, we discuss their thermal load models in details. We 
adapt a simple model, where the output temperature depends 
only on the current state k in which the device is operated i.e., 

hiiHii*)) = 

1) Thermal Model of AC: The output of an AC is the 
room temperature that can be obtained by its operation. The 
room temperature variation between any two consecutive time 
intervals is modeled by the following equation (inspired from 

USD, 

0*ac(* + !) ~0jUcW = + At^^-w jAC (t) (9) 

In this model, Gj(t) is the heat gain rate of the house of 
customer j, which depends on heat gain coefficients of the 
walls, windows, roof, solar radiation, people and air change 
rate of the AC, inside and outside temperature difference, etc., 
(■ Gj(t ) is independent of state k), a is the energy required for 
a unit degree rise in room temperature and Z* AC is the cooling 
capacity of the AC when it is operating in state k. Cooling 
capacity of an AC is specified by the manufactures in kW, 


BTU/hr or Tons. The cooling capacity is a function of the 
power state in which AC is operated and can be modeled as, 

z* AC = EER x q* ti (t) (10) 


where EER is the Energy Efficiency Ratio, which is typically 
defined as the ratio of cooling capacity given in BTU/hr to the 
power input in Watts. The higher the EER rating, the greater is 
the performance. The US national appliance standards dictate 
all AC loads to have a minimum value of EER > 8.0 1451. For 
example, an AC with a cooling capacity of 1 Ton (equivalent 
to 12000 BTU/hr or 3.516 kW) with an EER of 8.0 will 
consume 1.5 kW power. 

2 ) Thermal Model of WH: Hot water temperature obtained 
from operating the WH for one time slot in any given power 
state k can be modeled by the following equation (inspired 
from f[36l ). 


wh(^ ‘ 


1) = 


+ 


^j,WH (*)(Vtank “ fr{t).At) ^inlet- 




Mink 

Aank.(^, 


WH 


(t)-O a (t)) 


Rtimk 


^tank 

At 

ktank 


(ID 


In this model, water temperature in the next time slot t + 1, 
by operating WH in power state k (which consumes gjL(t) 
amount of power) depends on water temperature at the start of 
the time interval, temperature of the inlet water (0iniet)» current 
room temperature ( 0 a (t )), water flow rate ( fr(t )) during time 
slot t, volume of the tank (Vtank)» tank surface area (A t ank) and 
the heat resistance of the water tank (iZtank)- In this model, 
there are three terms and some necessary conversions might 
be required to make the units of all these terms consistent. 

3) Algorithm for Thermostat Loads: The home controller 
of user j receives the aggregated load profile denoted by 
Xj -1 from the grid. The objective of Algorithm [2| again is 
to reduce the peak load, while respecting constraints 0, 
® and © as laid out in the customer’s terms of engage¬ 
ment with the grid. There are two inconvenience dimensions 
(severity and duration) for each thermostat load. Thermostat 
loads are considered in a sequential order. Each thermostat 
load can be denied operation at its full rated power for 
a maximum of tf ir = At number of time slots in its 
demanded operation interval. Let b^ denote the set of time 
indexes during which customer j has demanded thermostat 
load i. For each thermostat load, the algorithm determines 
tf 11 local peaks in the interval by*. Let, t™ x = max^. . h 
and t™ x ~ n ,n = 1,..., tf 1 * + 1 represent the time index 
corresponding to n-th largest local peak. The time index vector 
representing the local peaks arranged in descending order is 
then denoted by t^f = [*£f, ..., A The 

objective is to operate the thermostat load i of customer j 
in its lowest possible state at time index t™ x , followed by 

Pff- 1 until -1 , without violating the inconvenience 

severity constraint anywhere in the demanded interval b^. 


The algorithm starts by switching OFF the thermostat load 
i of customer j at all the time indexes in vector t^ ak (i.e., 
c] : i(t) = 1, Vt G t^ ak ). Once the operational states are fixed at 
the desired time indexes, the algorithm determines a sequence 
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of output temperatures denoted by Oj^(t),\/t E bj^. If the 
inconvenience severity constraint is not violated anywhere 
in the demanded time interval, the algorithm will terminate. 
Otherwise, the algorithm will increase the operational state 

max—t dur — 1 

of the thermostat load i at the time index t- { 1 (i.e., 

c i,^ m r tf_1 ) = ° and = d and re - 

computes the output temperature sequence. This process is 
repeated until the inconvenience severity constraint is satisfied 
everywhere in demanded operational interval of the customer. 
This algorithm ensures that at the highest local peak points, 
the thermostat load is either switched OFF or it is operating 
in the lowest possible power state. Since least preference is 
given to the lowest local peak points, the thermostat load 
might operate at the full rated power at these points. Thus, for 
some customers, the algorithm will achieve the inconvenience 
duration constraint with strict inequality, in order to satisfy 
inconvenience severity constraint in all of its demanded time 
sloti. 


Algorithm 2 Algorithm for Thermostat Loads 

Initialize: h = xy _1 - ^2 ieJ s p jti , C jt i 0,Vi E 

T j 

i T 

1: for i E 4 do 

2: Determine and arrange in descending order the time 

indexes of tf m local peaks in the aggregated pro¬ 
file h in the time interval b ^ denoted by: t^ ak = 

[-/.max j.max-1 .max—tf* r -li 

I l j,i l j,i ' ' ' l j,i -I ^ 

3: Initialize: k = 0, t + a = li a , a = 1, cj^(t) = 

4: while 1 0j,i(t) — 0j : i\ > AOj^ for any t E b^ do 

5: Cj^t + a) = 0 

6: Increment: k = k + 1 

7: + 1 

8: Determine Oj^ using the thermal load model 

9: if k = K'i then 

10: k = 0, a = a + 1 

11: end if 

12 : end while 

13: p- ?i = {Qj,i ®Cj^)l K j (Note: © denotes element-wise 

matrix multiplication) 

14: Update the aggregated load profile: h = h + p^ • 

15: end for 


step does not result in significant increase in the complexity. 
For instance, if we assume that every customer has only two 
thermostat loads (AC and WH), then there are only two device 
orders: AC followed by WH or WH followed by AC, and 
Algorithm [2] will only be executed twice. 

TABLE II 

Appliance power rating and usage patterns 


Flexible 

Load 

Power Rating 
(kW) 

Operation 

Duration 

Usage 

Pattern 

AC 

5 

4 hours 

Consecutive 

operation 

WH 

2.5 

3-6 hours 

Two or three separate 
instances of 1-2 hours 

CD 

3.1 

2 hours 

Consecutive 

operation 

DW 

1.8 

1.5 hours 

Consecutive 

operation 


IV. Simulation Results 

We consider a residential community comprising of 1000 
homes. Each customer is assumed to have two thermostat loads 
i.e., AC and WH and two shiftable loads i.e., CD and DW. All 
other devices contribute towards an essential base load. The 
average daily household energy consumption is assumed to be 
about 41 kWh, which corresponds to typical household energy 
consumption in many US states like Louisiana, Tennessee, 
Alabama etc. The appliance usage and the power consumption 
pattern of each device is given in Table II, which are generated 
using realistic assumptions as given in 106), 03 and also to 
match a load curve shape from the RELOAD database on a 
typical summer day. The RELOAD database is an industry 
accepted database of load curve shapes and is used by the 
electricity module of the National Energy Modeling System 
(NEMS) and several authors for their studies ||48l, [|49lL ll36l . 
The thermostat set points of AC and WH are uniformly 
distributed random variables in the interval [68°F,76°F] and 
[104°F,120°F]. The resulting aggregated load profile of the 
community is shown in Fig. [3] with a peak load value of 3701 
kW. In the simulations, we assume Kj = K , Vj, i E I° T and 
evaluate the impact of 2-state, 3-state, and 5-state models for 
AC and WH. The EER value of AC load is assumed to be 
10. We run the simulations in MATLAB and divide the total 
time duration i.e., 24 hours into T -288 equal interval time 
segments (A t = 5 mins). 


This algorithm again requires optimization over all devices 
orders, since considering thermostat loads in a different se¬ 
quence can result in a different amount of peak load reduction. 
This can be done by executing Algorithm [2] for (I J T )\ times 
at home controller of customer j. The number of thermostat 
loads in every home is limited, therefore, this optimization 

5 The worst case complexity in terms of number of iterations is 
ieI j (7C)!. For example, if we consider two thermostat loads per customer 
each having 5-states, the worst case complexity of Algorithm [2] is 240 
iterations. We show in the simulations, that increasing the power states beyond 
three is not much beneficial for peak load reduction. The number of thermostat 
loads as well as the power throttling states of devices are generally low, 
therefore, this complexity is manageable. 



Time (Hours) 

Fig. 3. Aggregate load profile of the residential community 

We evaluate the effectiveness of CDPs and PDPs and the 
impact of increasing the power throttling states on peak load 
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Fig. 4. Percentage peak reduction achieved by CDPs. Y-axis shows the 
maximum inconvenience durations of AC and WH loads in the format: 

A^ac* -The maximum temperature deviations of AC and WH 

loads are shown along the bars in the format: (AO™^ -A6^[). 
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Fig. 5. Percentage peak reduction achieved by PDPs. Y-axis shows the 
maximum inconvenience durations of AC and WH loads in the format: 

A^a? -Ad^[. The scaling factors of AC and WH loads are shown along 

the bars in the format: (J3ac - Avh)- 


reduction in Figs. 0] and \5\ respectively. The common plan 
parameters in both these figures are: Asj^ x = 0 mins and 
A = 0 mins, Of c = 80°F, 6™ f n = 96°F. In Fig. ffl we 
vary the parameters, which define the temperature deviations 
and inconvenience durations of AC and WH loads for CDPs. 
While in Fig. 0 we vary the scaling factors and inconvenience 
durations of AC and WH loads for PDPs. Based on these two 
figures, we have the following observations. 

• Increasing the inconvenience durations of AC and WH 
loads, the temperature deviations in CDPs, the scaling 
factors in PDPs, and the number of states of AC and 
WH loads increase the peak load reduction. 

• When inconvenience durations of AC and WH loads in 
CDPs and PDPs are fixed, we observe a marginal peak 
load reduction with increasing temperature deviations 
in CDPs and scaling factors in PDPs (also termed as 
diminishing returns). This is due to the saturation of the 
inconvenience severity dimension of the plan. Similarly, 
when temperature deviations of AC and WH loads in 
CDPs and scaling factors in PDPs are fixed, we again 
observe diminishing returns when the inconvenience du¬ 
rations of loads are increased. This is due to the saturation 
of the inconvenience duration dimension of the plan. 

• The inconvenience severity dimension saturates at a faster 
rate compared to the inconvenience duration dimension 
of the plan for both the CDPs and PDPs. 

• For both CDPs and PDPs, diminishing returns are ob¬ 
served with an increase in the number of power throttling 
states. 


• There is a tradeoff between the returns on the number 
of states and inconvenience severity dimension of the 
plan. Plans offering less inconvenience severity (less 
temperature deviations for CDPs and less scaling factors 
for PDPs) exhibit more gains with the increase in the 
number of power states. On the other hand, plans with 
high inconvenience severity, do not offer any significant 
gains with the increase in the number of states. In 
high inconvenience severity plans, intermediate states 
are generally not required even if they are available, 
since operating the loads in OFF states save more power 
without violating the inconvenience severity constraints. 

In Fig. 0 we plot the impact of scheduling delays of CD 
and DW loads on the peak load reduction. The CDP plan 
parameters are: AC load (60 min, 2°F, 80°F), WH load (60 
min, 4°F, 96°F), while PDP parameters are: AC load (60 min, 
0.25, 80°F), WH load (60 min, 0.25, 96°F). The temperature 
deviations experienced by the AC and WH loads of different 
customers are different in PDP but constant in CDP, which 
can be visualized from the top graph in this figure. In these 
simulations, the choice of /3ac = Avh = 0.25 in PDPs, results 
in an average temperature deviation of 2.04°F for AC load and 
4°F for WH load, which match the temperature deviations of 
2°F for AC and 4°F for WH in CDPs. From the bottom graphs, 
we can observe that as we increase the scheduling delays, 
the amount of peak load reduction and the rate of peak load 
reduction both increase. Furthermore, PDPs, while offering 
more fair temperature deviations (inversely proportional to 
thermostat set points), also achieve almost identical peak load 
reduction that is obtained by the CDPs. 


Temperature Deviations of AC and WH loads in CDP and PDP) 



Temperature Deviation of AC in PDP 
“ “ “ Temperature Deviation of WH in PDP 
- ■- Temperature Deviation of AC in CDP 
~ 0~ Temperature Deviation of WH in CDP 



Scheduling Delays of CD and DW Loads in minutes 


Fig. 6. The top graph shows the temperature deviations experienced by the 
AC and WH loads of 100 customers in the community in CDP and PDP. The 
bottom graphs show the percentage peak reduction for different scheduling 
delays for CD and DW loads in CDPs and PDPs. X-axis of the bottom graph 
shows the maximum scheduling delays of CD and DW loads in minutes, in 
the format: As™g x -As^JJ. 


Finally, we study the impact of reference temperature on 
peak reduction for CDPs and PDPs and the results are pre¬ 
sented in Table III. Common Parameters of CDPs and PDPs 
are As™ x = A= A d™ x = Ad™* = 60 min. CDP 
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TABLE III 

This table shows the impact of reference temperatures on the number of eligible loads, average temperature deviations over 

THE DEMANDED INTERVAL AND % PEAK REDUCTION FOR CDPS AND PDPS. 


Ref. Temp. 

No. of eligible loads 

Average Temperature Deviations 

% Peak Reduction 

mve 

'Ac 

mve 

P WH 

2-state 

3-state 

5-state 


nref 

P WH 

Aac 

iV WH 

CDP 

PDP 

CDP 

PDP 

CDP 

PDP 

CDP 

PDP 

CDP 

PDP 

72°F 

112°F 

550 

509 

1.51°F 

0.5°F 

2.97°F 

1.0°F 

14.98 

6.13 

22.66 

20.71 

24.49 

23.45 

74°F 

108°F 

766 

752 

1.71°F 

0.79°F 

3.32°F 

1.50°F 

19.36 

8.07 

23.95 

21.27 

25.4 

24.28 

76°F 

104°F 

1000 

1000 

1.76°F 

1.04 

3.48 

2.0 

22.84 

13.55 

25.45 

22.58 

26.22 

25.03 

78°F 

100°F 

1000 

1000 

2°F 

1.54°F 

4°F 

3°F 

24.1 

21.22 

26.09 

24.66 

26.52 

25.91 

80°F 

96°F 

1000 

1000 

2°F 

2.04°F 

4°F 

4°F 

24.1 

23.78 

26.09 

26.08 

26.52 

26.51 


parameters: A/9™£ x = 2°F, A6 ^h = 4°F. PDP parameters: 
Pac = Avh = 0.25. In this table, N A c and Awh respectively 
denotes the number of AC and WH loads with non-zero in¬ 
convenience severity, while and respectively denotes 
the average temperature deviations experienced by Nac and 
Awh loads. Different values of reference temperatures result in 
different values of Nac and Awh- As we increase the reference 
temperature for AC and decrease that of WH, more peak load 
reduction is observed for CDPs and PDPs, since more devices 
are controlled. We can observe another tradeoff between the 
number of states and reference temperatures (which impact 
Aac and Awh) especially for PDPs. When the reference tem¬ 
peratures result in less number of ACs and WHs experiencing 
inconvenience, providing an additional power throttling state 
can provide very high performance gains compared to the 2- 
state (basic ON/OFF switching) models. We can also observe a 
similar trend (but in lesser magnitude) for CDPs. Furthermore, 
when 0$, > § jtAC + AVj and < % WH - Al9$g,Vj, 
then increasing 0 and decreasing 6 >Uh has no further impact 
on the peak load reduction in CDPs (since there is no further 
impact on the number and temperature deviations, as evident 
from the last two rows of Table III). On the other hand, in 
PDPs, increasing 0 r ^ c and decreasing 6^ f H always increase the 
amount of peak load reduction (since temperature deviations 
also increase). 

V. Conclusions and Future Research Directions 

In this paper, we have designed two types of customer 
engagement plan, namely CDP and PDP, which describe the 
key inconvenience parameters of flexible loads in terms of 
scheduling delays and temperature deviations so as to make 
the customers easily understand the inconvenience caused by 
these plans. To facilitate the grid operator in determining 
the effectiveness of such plans on peak load reduction, we 
have developed appropriate DRP algorithms by modeling 
the thermostat loads as power throttling devices. Despite the 
suboptimality (due to NP nature of the optimization problem), 
the proposed algorithms have been shown to be able to provide 
very clear insights into the design of customer engagement 
plans. Through simulations, we have determined that increas¬ 
ing temperature deviations result in diminishing returns. We 
have also observed diminishing returns as the number of 
power throttling states increase. In particular, more peak load 


reduction occurs when the number of power throttling states 
are increased from 2 to 3 for those DRP plans with less 
number of thermostat loads. The temperature deviation from 
thermostat set point cannot exceed the given limit in all the 
time slots, therefore, with two states, the only option to avoid 
exceeding the temperature deviation limit is to turn ON the 
appliance (consuming full power). On the other hand, the third 
state provides more flexibility and we can avoid exceeding the 
temperature deviation limit by operating the device at 50% 
power (hence resulting in 50% power reduction) compared to 
the two state model. 

We also have some recommendations for the design of cus¬ 
tomer engagement plans. The grid operator should generally 
design plans with low to moderate inconvenience severity. In 
a community with thermostat loads having more fine control 
(three or more states), the plans that offer low inconvenience 
severities are more beneficial. On the other hand, for a 
residential community with only ON/OFF power control for 
thermostat loads, more peak reduction can be achieved by de¬ 
signing plans with high inconvenience duration of thermostat 
loads and scheduling delay of shiftable loads. PDPs can be 
offered to the customers of a residential with highly variable 
thermostat set point preferences, while CDPs can be offered 
to the customers with similar thermostat set point preferences. 
The grid operator can use the reference temperatures in the 
plans to control the number of eligible devices. The value 
of reference temperatures can also be adjusted by the grid 
operator depending on the weather patterns, geographical 
location, user preferences, surveys etc. 

These results can be seen as useful guidelines in proposing 
financial rewards and mutually beneficial incentives, which 
is an interesting future work. Extending the framework to 
include more load categories e.g., electrical vehicles are an¬ 
other interesting research direction. Modeling the transient 
behavior and reducing the additional power consumption due 
to frequent switching of thermostat appliances from one power 
throttling state to another might also be an interesting future 
work. Finally, grid operator may also introduce multiple plans, 
each with different parameters, to cater for the preference of 
different customers. Managing and designing various plans 
also appears to be an interesting future work. 
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